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Abstract. We derive an expression for the entropy of a present dark matter halo described by a 
Navarro-Frenk- White modified model with a central core. The comparison of this entropy with the 
one of the halo at the freeze-out era allows us to obtain an expression for the relic abundance of 
neutralinos, which in turn is used to constrain the parameter space in mSUGRA models, when used 
with the WMAP observations. Moreover, by joning these results with the ones obtained from the 
usual abundance criteria, we are able to clearly discriminate validity regions among tan /3 values 
of the mSUGRA model, by demanding both criteria to be consistent with the 2 sigma bounds of 
the WMAP observations for the rehc density: 0.112 < f^h'^ < 0.122. We found that for sgn p = +, 
small values of tan/3 are not favored; only for tan/3 ~ 50 are both criteria significantly consistent. 
The use of both criteria also allows us to put a lower bound on the neutralino mass, > 151GeV. 

PACS. 14.8Q.Ly Supersymmetric partners of known particles - 95.35.+d Dark matter - 98.62.Gq 
Galactic halos 



1 Introduction 

Supersymmetry models which have the neutrahno as 
the lightest supersymmetric particle (LSP) and as a 
candidate for dark matter (DM), have several param- 
eters that can be constrained by the bounds on the 
present density of DM, ficDM, that come from sev- 
eral outstanding observations such as the Cosmic Mi- 
crowave Background radiation (CMBR) [1], Galaxy 
clustering, Supernovae and Lyman a forest. One of 
the most recent works which combines all these data 
leads to: 0.112 < QcDKih^ < 0.122 [2 . In particu- 
lar for mSUGRA models these constraints have been 
obtained using the standard approach [Sll^, which is 
based in the Boltzmann equation considering that af- 
ter the "freeze-out" era, neutralinos cease to annihi- 
late keeping its number constant. In such an approach, 
the relic density of neutralinos is approximately: « 
l/((Tf), where (av) is the thermally averaged cross sec- 
tion times the relative velocity of the LSP annihilation 
pair. Within the mSUGRA model five parameters (toqi 
'7^-1/2, ^Oi tan/3 and the sign of fi) are needed to spec- 
ify the supersymmetric spectrum of particles and the 
final relic density. We will use the numerical code mi- 
crOMEGAs to compute the relic density following 
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the past scheme which will be called the "abundance 
criterion" (AC). 

Just after "freeze-out" , we can consider neutrali- 
nos then as forming a MaxwcU-Boltzmann (MB) gas 
in thermal equilibrium with other components of the 
primordial cosmic structures. In the present time, such 
a gas is almost colisionless and either constitutes galac- 
tic halos and larger structures or it is in the process 
of its formation. In this context, we can conceive two 
equilibrium states for the neutralino gas, the decou- 
pling (or "freeze-out") epoch and its present state as 
a virialized system. Computing the entropy per parti- 
cle for each one of this states we can use an "entropy 
consistency" criterion (EC) using theoretical and em- 
pirical estimates for this entropy to obtain the relic 
density of neutralinos (i^^)- 

Our objective is then to use AC and EC criteria, to 
obtain constraints for the parameters of the mSUGRA 
model by demanding that both criteria must be con- 
sistent with each other and with the observational con- 
straints of ficDM- 



2 Abundance criterion 

Rehc abundance of some stable SUSY species x is de- 
fined as f2^ = Pxl Pcrit, where = m^n^ is the relic's 
mass density [n^ is the number density), Pcrit is the 
critical density of the Universe (see [HI for a review 



Cosmology and Astrophysics 



Contributed Talk 



on the standard method to compute the reHc density) . 
The time evolution of is given by the Boltzmann 
equation: 



dn^ 



{av){n 



(1) 



where H is the Hubble expansion rate, (av) is the ther- 
mally averaged cross section times the relative velocity 
of the LSP annihilation pair and n^' is the number 
density that species would have in thermal equilib- 
rium. In the early Universe, the neutralinos (x) were 
initially in thermal equilibrium, ~ n'^ . As the Uni- 
verse expanded, their typical interaction rate started 
to diminish an the process of annihilation froze out. 
Since then, the number of neutralinos in a comoving 
volume has remained basically constant. 

There are several ways to solve equation ([T]), one 
of the more used is based on the "freeze-out" approxi- 
mation (see for example j^)- However in order to have 
more precision, we will use the exact solution to Boltz- 
mann equation using the public numerical code mi- 
crOMEGAs 1.3.6 [5 which calculates the relic density 
of the LSP in the Minimal Supersymmetric Standard 
Model (MSSM). We will take and mSUGRA model 
and its five parametersas input parameters for mi- 
crOMEGAs and use Suspect [H^, which comes as an 
interface to micrOMEGAs, to calculate the supersym- 
metric mass spectrum. 

Using micrOMEGAs, we can obtain the relic den- 
sity for any region of the parameter space to discrimi- 
nate regions that are consistent with the WMAP con- 
straints in this abundance criterion. 



3 Entropy consistency criterion 

Since the usual MB statistics that can be formally ap- 
plied to the neutralino gas at the "freeze-out" era can 
not be used to describe present day neutralinos sub- 
ject to a long range gravitational interaction making 
up non-extensive systems, it is necessary to use the ap- 
propriate approach that follows from the microcanoni- 
cal ensemble in the "mean field" approximation, which 
yields an entropy definition that is well defined for 
a self-gravitating gas in an intermediate state. Such 
an approach is valid at both the initial ("freeze-out" 
era, /) and final (virialized halo structures, h) states 
that we wish to compare. Under these conditions, the 
change in the entropy per particle (s) between these 
two states is given by j9]: 



In 



rf 



3/2' 



(2) 



where x — m^/T, T is the temperature of the gas. 
A region that fits with the conditions associated with 
the intermediate scale is the central region of halos 
( IQpc^ within the halo core); evaluating the thermo- 
dynamical quantities at this region, using equation ^ 
and some extra assumptions (conservation of photon 



entropy), it is possible to construct a theoretical esti- 
mate for that depends on the nature of neutralinos 
(m^ and {ov)), initial conditions (given by x^), cos- 
mological parameters (J^^, the Hubble parameter, h) 
and structural parameters of the virialized halo (cen- 
tral values for temperature and density); for details of 
these and the following, see section IV of [5]. 

An alternative estimate for s'' can be made based 
on empirical quantities for observed structures in the 
present Universe using the microcanonical entropy def- 
inition in terms of phase space volume, but restricting 
this volume to the actual range of velocities accessi- 
ble to the central particles. That is, restricting the es- 
cape velocity up to a maximal value We(0) which is 
related to the central velocity dispersion of the halo 
(cT/i) by an intrinsic parameter a: v1{Q) ~ acr^(O). In 
a recent work [10', we estimate the value of a using 
an NEW modified model with a central core, and ob- 
tain 16.4 < a < 27.8. The range of values allowed for 
this parameter is of the highest importance to deter- 
mine the allowed region of the parameter space in the 
mSUGRA model as will be clear in the results pre- 
sented on next section. 

Equating the theoretical an empirical estimates for 
the entropy per particle it is obtained a relation for the 
relic abundance of neutralinos using the EC criteriorj^ 



ln{n^h^) 
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where fg{x-f) is a function related to the degrees of 
freedom at the "freeze-out" time (see for example [7]) 
that will be described elsewhere [TU . 

Modifying the program micrOMEGAs, we obtained 
the value for x^ for any region of the parameter space 
and then Q-^ using Eq. ([3]), therefore we were able 
to discriminate regions that are consistent with the 
WMAP constraints for the EC criterion. 



4 Results and Conclusions 

Using both the AC and EC that have been described in 
the preceeding sections, we can compute the total mass 
density of neutralinos present today and constrain the 
region in the mSUGRA parameter space where both 
criteria are fulfilled. Out of the five parameters, we will 
fix /i > 0. Then our strategy is to explore wide regions 
for the values of the other four parameters, by means 
of a bi-dimensional analysis in the mp — mi/2 plane for 
different fixed values of Aq and tan/3. It is important 
to mention that we are not presenting an exhaustive 
search in all the possible regions, but we concentrated 
on those regions which have received more attention 
in the literature, see for example [3]. 

In Fig. ([ij, we present the results for tan/3 = 10, 
for three values of ^o, namely Aq — 1000,0,-1000 
GeV, shown in the top, middle and bottom panels re- 
spectiveley. The yellow region (lower right corner) is 



^ This formula is a small modification to the one pre- 
sented in [9] 
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Fig. 1. Allowed regions in the parameter space for AC 
(lighter gray /red) and EC (darker grey/blue) for the 
mSUGRA model with sgn/i = +, tan/3 = 10, and Aq — 
1000 GeV, top panel, Ao — GeV, middle panel, and 
Ao = —1000 GeV, bottom panel. The figures show the so 
called bulk and coannihilation regions. The yellow region 
shows where the stau is the LSP. 



where the f is the LSP, the lighter and darker areas 
(red and blue for the online version in colours) define 
the allowed regions for the EC and AC respectively ac- 
cording to the observed DM density. The area of the 
EC region depends on the size of the interval of val- 
ues of the parameter a, the lower and upper bounds 
of a determine the upper and lower boundaries of the 
EC region. As can be seen from the figure, the region 
where both criteria are fuUfiUed is very small, in fact, 
only for the highest values of a there is an intersection 
between both criteria. This behavior holds for all val- 
ues of Aq in the interval [-1000,1000] GeV, here we 
are showing only the extreme and central values. 

Repeating the same procedure for larger values of 
tan/3, it is found that the intersection region for both 
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Fig. 2. The same as Fig. |l]), but for tan/3 = 50, and now 
Aq — —500 GeV in the bottom panel. 



criteria becomes larger, but it gets to be significant 
for the largest values of this parameter. This is clearly 
shown in Fig. which is equivalent to Fig. (jlj, but 
for tan/3 = 50. In this case the bottom panel is for 
^0 = —500 GeV. It is clear from the figure that for 
these values of tan/3 both criteria are consistent, as 
shown by the large intersection area for values of Aq 
in the interval [0, 1000] GeV. For negative values of 
Aq the intersection region decreases with Aq, see the 
bottom panel of the figure. For even lower values of Aq 
the intersection becomes insignificant. 

In Figs. ([3]) we present the same analysis but for the 
Focus Point region, and for the central value ^0 — 0. 
The situation is consistent with the previous results, 
both criteria intersect for tan /3 = 50 and there is 
nearly no intersection for tan/3 — 10. 

This analysis allows us to arrive to one of the main 
results of our work. The use of both criteria favours 
large values of tan/3. 
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Fig. 3. Allowed regions in the parameter space for 
AC (lighter gray/red) and EC (darker grey/blue) in the 
mSUGRA model with Aq — 0, sgn/i — +, tan/3 — 10, top 
panel, and tan/3 = 50, bottom panel. The region shows the 
so called Focus Point region. 



In Figs.(|4]) and ([5]) we show the allowed values for 
the LSP ana the Higgs mass after constraining the 
parameter space with the abundance and entropy cri- 
teria. As can be seen from Fig. (|4]), the current limit 
for the Higgs favour, combined with the AC and EC 
criteria, favours even more a large value of tan /3. This, 
in turn, puts a constraint on the allowed SUSY mass 
spectra of the bulk and coannihilation regions: it gives 
an LSP of mass ^ 140 GeV for tan (3 10, and a 
lower bound for the LSP mass > 150 GeV for 
large tan/3. 

Further analysis, which is currently under way, is 
required to give more precise conclusions about this 
new method to constrain the parameter space of the 
mSUGRA model |l0]. 
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Fig. 4. Allowed values for Muiggs as function of m-i/2- 
As can be seen from the figure, the present bound on the 
Higgs mass favours a large tan /3. 
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Fig. 5. The figure shows the LSP mass plotted versus the 
Higgs mass, points above the dashed line are the allowed 
values for the LSP. The points in blue correspond to the 
Focus Point region, the ones in red to the bulk and coan- 
nihilation regions. 
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